Radiation back-reaction in relativistically strong and QED-strong laser fields 
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The emission from an electron in the field of a relativistically strong laser pulse is analyzed. At 
the pulse intensities of > 10 22 W/cm 2 the emission from counter-propagating electrons is modified 
by the effects of Quantum ElectroDynamics (QED), as long as the electron energy is sufficiently 
high: £ > 1 GeV. The radiation force experienced by an electron is for the first time derived from 
the QED principles and its applicability range is extended towards the QED-strong fields. 



In Quantum ElectroDynamics (QED) an electric field 
should be treated as strong if it exceeds the Schwinger 
limit: Eg = m e c 2 / (\e\Xc) PJ. Such field is potentially ca- 
pable of separating a virtual electron-positron pair pro- 
viding an energy which exceeds the electron rest mass 
energy, m e c 2 , to a charge, e, over an acceleration length 
as small as the Compton wavelength, Ac = • Partic- 
ularly a QED-strong Lorentz-transformed electric field, 
Eo = |px B|/(m e c), may be exerted by a charged particle 
with the momentum p, gyrating in the magnetic field, B, 
if |p| 3> m e c and/or the magnetic field is strong enough. 

Consider QED-effects in a strong pulsed laser field [3] : 



x a 2 > 1, 



eA 



m e c 



(1) 



A being the vector potential of the wave. In the labo- 
ratory frame of reference the electric field is not QED 



3], 
4]. 



strong for achieved laser intensities, ~ 10 22 W/cm 2 
and even for the ~ 10 25 W/cm 2 intensity projected 
Nonetheless, a counter-propagating particle in a ID wave, 
a(£), £ = u>t— (k-x), may experience a QED-strong field, 
Eo = \dA/d£\uj(£ —pn)/c, because the laser frequency, 
u), is Doppler upshifted in the frame of reference comov- 
ing with the electron. Herewith the electron dimension- 
less energy, £, and its momentum are related to m e c 2 , 
and m e c correspondingly and subscript || herewith de- 
notes the vector projection on the direction of the wave 
propagation. The Lorentz-transformed field exceeds the 



Eo/Es 



> i, 



Schwinger limit, if \ 
where X = — . 

to 

Within classical theory a signature for Eq is a radia- 
tion loss rate, I C \{E ) = 2e 4 E 2 / (3m 2 c 3 ). Therefore, the 



QED-strcngth of the electromagnetic field may be deter- 
mined in evaluating I c \ and its ratio to Ic — I c \{2Es /S): 




Ic = 



27X7 



(2) 



'Electronic address: igorsok@umich.edu 



If X ^ 1 then the actual radiation loss rate differs from 
7 c l. The condition of \ > 1 also separates the parameter 
range of the Compton effect from that of the Thomson 
effect, under the condition of Eq.|T|). The distinctive fea- 
ture of the Compton effect is an electron recoil, which 
is significant, if a typical emitted photon energy, Huj c , 
is comparable with the electron energy ||. Their ratio, 
X = XcuJcl {c£), equals x as defined in Eq.((2]) with the 
proper numerical factors (cf Eq. p^]) ). 

Counter-propagating electrons can be generated in the 
course of laser pulse interaction with a solid target, that 
is why the radiation effects in the course of laser-plasma 
interaction are widely investigated (see @, @). The prin- 
ciple matter in this paper is an account of the radia- 
tion back-reaction acting on a charged particle. This 
can be consistently done by solving the modified Lorentz- 
Abraham-Dirac equation as derived in [7J, in which the 
radiation back-reaction on the electron motion is ex- 
pressed in terms of the emission probability. The cal- 
culation of this probability in relativistically strong and 
QED-strong laser pulses is given in Section I. In Section 
II we discuss the radiation effect on the electron motion 
in strong fields. 



I. ELECTRON IN A ID WAVE: THE EMISSION 
PROBABILITY 

The emission probability in the strong ID wave field 
may be found in §§40,90,101 in 0, as well as in 
However, to simulate highly dynamical effects in pulsed 
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FIG. 1: The volume over which to integrate the matrix ele- 
ment while finding the emission probability: in the standard 
scheme for the dipole emission (in dashed lines) and in the 
TST (in solid lines). Arrows show the direction, along which 
the Heisenberg operator advances the wave functions. 



fields, one needs a reformulated emission probability, re- 
lated to short time intervals (not (—00, +00)). 

Consider the classical motion of an electron in 

a ID field, a = a M (£), £ = (k ■ x), (k ■ a) — 0, a, k and 
x being the 4-vectors of the potential, the wave and the 
coordinates. Herewith the 4-dot-product is introduced in 
a usual manner: (k ■ x) = Lot — (k • x) etc., 3- vectors in 
contrast with 4-vectors being denoted in bold, 4-indices 
are omitted below. Introduce a Transformed Space-Time 
(TST) : x ' 1 = (ct T x\\)/V2, x 2 ' 3 = x±, subscript _L 
denoting the vector components orthogonal to k. Note, 
that: (1) dx° = Xd£/V2, p° = X{k- p)/V2, (p-k) = (E- 
P\\)/X; (2) the momentum components, p° and + a, 
are conserved; (3) the metric tensor in the TST is: G 01 = 
G 10 = 1, G 22 = G 33 = -1; and (4) p^G^p" = 1 gives: 
p 1 = (1 + p 2 1 )/(2p°). These properties allow us to find: 

2e 2 cdpj rjk dp k _ 2e 2 c(k-p) 2 ' ' * 1 
ds 



Id = 



da 



3 ds ds 3 
and to relate 4-momenta at different time instants 



(3) 



K0=p(r) _ {o+ wa_M^ t , (4) 



where ds = y / c 2 dt 2 — (p • dx/£) 2 and 5a = a(£) — a(£'). 

A QED solution of the Dirac equation in the TST 

is given by a plane electron wave (N = const), 



*(p(0) 



N 



C(£) exp 



«(pio-xi) i(k-p)x l 



Xc 



(5) 



(A C /A)V2 

The Dirac equation, ((7 • (iXc-§^ — a)) — l) i/) = 0, is sat- 
isfied under the following conditions: u(p(£)) is plane 
wave bi-spinor amplitude, ((7 • p(£)) — 1) u(p(£)) = 0, 
u ■ u — 2, 7^ are the Dirac matrices, and C(£) = 



exp i 



«(p(0) = 



(7 ■ k) (7 ■ 5a) 
2{k-p) 



(6) 



The emission probability. Introduce domain, 
A 4 a; = (A3; 1 * S±) * Ax° , bounded by two hypersur- 
faces, £ = £- an( i £ = £+ (see Fig.l). The differ- 
ence £+ — £_ is bounded as decribed below, so that 
A 4 x covers only a minor part of the pulse. A volume 
V = S±X(£ + — £_), is a section of A 4 x subtended by a line 
t = const. With the choice of the coefficient in Eq.©, 

N = 2S±X £(£ 2 )d£ 2 , the integral S±X /|+ V^V^i 
is set to unity, i.e. there is a single electron in the vol- 
ume V. The emission probability, dW, for a photon of 
wave vector, k' , polarization vector, I, and wave function, 



A' 



i{k' -x)] 



u'V 
2irhc' 2 ■ 



is given by an integral 



over A 4 a;: 



dW 



uL fLp 
he 



tp f (j ■ (A')*)%l} i dx dx 1 dx 2 dx i 



(J) 



L p = d?k' x d(k ■ k')hcN p /((2Tr) 2 (k ■ k')) is the number of 
states for the emitted photon, a subscript i, / denotes the 
electron in the intitial (i) or final (f ) state, and a = f- « 
j^j. The number of electron states, Lij, per d 3 Pij = 
d 2 p±ijd(k ■ p)ij£ij (£) / (k • p)ij , in the wave field can be 
found by calculating the Hamiltonian invariant, d 3 pdV ~ 

g(£)d£, with no field: L hf = M^^ , ^ 

Conservation laws. The integration by da; 1 da; 2 da; 3 = 
c-*/2dtd 2 x± results in three 5— functions, expressing the 
conservation of totals of and (k ■ p), for particles in 
initial and final states. Twice integrated with respect to 
dx 1 , the probability dW is proportional to a long time 
interval, At = Ax 1 / (c\/2), if the boundary condition for 
the electron wave at £ = £_ is maintained within that 
long time. For a single electron, which locates between 
the wave fronts £ = £_ and £ = £ + during a shorter time, 



«ft(£_,£+) = (l/c) / £i(Odea/(fc-Pi), 



(8) 



the emission probability is: dWfi (£_,£+) = <5tdVF/At. 
Using <5— functions we integrate Eq.([7|) over dp±fd(k-pf): 



d(fc • fc')d 2 k' ± ~ (47r) 2 (fc • k'){k-pi)(k-p f ) 
where F(£) = exp[z J ? (*' • Pifa))d&/(k -p f )] and 

(*'-Pi(0)*c 



(k ■ pi) - (k ■ k')X c 



k. 



(9) 



(10) 



To integrate Eq.([9]), we re- write it as the double 
integral over ei£d£i and reduce the matrices u(pij(£)) <g> 
^(Pi,/(Ci)) m trie integrand to the polarization matrices 
of the electron at £ or at £1 using Eq.©. Although in 
a strong wave electrons may be polarized (see [Til]), we 
then average over electron and photon polarizations and 
find: 



dW 



fi 



* it ft mn-ti)Mz,zi)d£dti 



d{k ■ k')d 2 k' ± 



{2Tr) 2 {k-k')(k-pi)(k- Pf ) 
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where D = 



1 . Now we 



^)-a(C 1 )] 2 ((fc.p,) 2 + (fc.p / ) 2 
i(k- Pi )(k-p f ) 

develop the dot-product, (k'-pi), in F(£) in the TST met- 
ric G and find: F(£)F(-£i) = exp [i{F x + F 2 )], where 



Ft 



(k ■ Pi 



(k ■ k') 

2(k-k')(k- Pf ) \Tk~p~) {p±l 

(fc-fc') {(?-?!)+ [a(C 2 )-(a)] 2 d5 2 | 
2(fc- Pl )(fe-p/) 



f 5 ad? 2 



and 



d(k ■ k') 



2ir(k -pi) 2 



In the strong field as in Eq.([T]) the formulae sim- 
plify. In F2 we estimate: £ — £1 ~ |da/d£| , and 
(fc • k') ~ (fc • Pi)(fc • Pi) |da/d£|. Now we can conistently 
introduce the bounds for £+ — £-: 



|da/d£| 1 <£+ -£_ <min (aT 1 |da/d£| \l 



Under these bounds, the emission probability: (1) is lin- 
car in £+-£_: dW fi (£_,£+) = (dW/i/d£)(£+ - £_); 
(2) is less than unity: / dW/j < 1; and (3) can be 
expressed in terms of the local electric field. By in- 
traducing (0 = (£ + £0/2, 9 = (£ - £ 1 )|da/d£|/2, 
(a(£) — a(£i)) 2 ss 46* 2 and expressing the integral over 
9 in terms of the MacDonald functions wc find: 



dW, 



■ fi _ a X (/" #f (y)dy + r r xX 2 K2 (r x )) 



drod£ 



%/37r(fc • k)Ac 



(12) 



X = -^{k-pi) 



da 

d£ 



Id 
Ic 



1 - xro ' 

Probability (similar to that found in 0) is expressed in 



terms of functions of ro = 



Qfc') 



and related to 



3 (k-pi) 2 \da/d£\ ' 

interval of dro . Below we demonstrate the way to use this 
probability to describe the electron motion and emission. 



II. RADIATION AND ITS BACK-REACTION 

At 3> 1 the wave vector of the emitted photon is 
almost parallel to the electron momentum. For colinear 
k' and ft, one has oj'/£i « c(k ■ k')/(k ■ p%), therefore, 



r Q = — , uj c 



-Sc(k-pi) 



da 

d£ 



£c X /X 



c- 



(13) 



The assumption, 3> 1 also allows us to find the mo- 
mentum of the emitted radiation, which we relate to the 
interval of the electron proper time, using Eqs. (|8ll2[) : 



dp T 



ad 



dr 



= / Xr:k' 



c(k ■ Pi )dW 
d{k • fc')d 2 k ± d£ 



d(k ■ k')d 2 k ± = 
dW 



drod£ 



<y> —4 
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FIG. 2: Emission spectra for various values of \- 
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FIG. 3: Emitted radiation power in the QED approach vs 
classical (solid); an interpolation formula Tqed = J c i/(1 + 
l.OV/d/Ic) 473 (dashed). 



As with other 4-momenta, p ra <j is related to m e c. To 
prove the 4- vector relationship (|14|) . we expand its com- 
ponents in the TST metric and integrate them over kj_ 
using the symmetry of F\. The small term, 0(l/(kpi)), 
arises from the electron rest mass energy and from the 
small (~ l/|pj_|) but finite width of the photon an- 
gular distribution. Below we neglect this term and 
find: dp rad /dr = Pi^-^(dlQED/dr )dr where / QE d = 



dW 
d^dro 



drn is the radiation loss rate. 



m e c 2 J cXc(kk') 

Thus, the angular distribution can be represented as 
5(fl — p/ \J~p 2 ~)dQ., with Q being the solid angle of 
the photon direction. The photon energy spectrum, 
d/QED/dro, is described as function only of the random 
scalar, ro = to' /cj c , using only the parameter, \ ( see 
Fig©- The latter may be parameterized in terms of the 
radiation loss rate, evaluated within the framework of 
classical theory (see Eq.([2]) and Fig(3j). The expressions 
for q{I c \) = Iqet>/Ic\ and for the normalized spectrum 
function, Q(tq,x)i coincide with formulae known from 
the gyrosynchrotron emission theory (see §90 in [8j]): 



9^ 



[p + k 0((k ■ p^- 1 )] / cXc(k'-k)— ^—dr . (14) Q(r ,x) 



dr r 



9\/3> 
8irq 



Ks. {y)dy + r r x x 2 K2 (r x ) 



K5 (y)dy + r r x x 2 K % (r x ) 
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FIG. 4: The emission spectrum for 600 MeV electrons inter- 
acting with 30-fs laser pulses of intensity 2 • 10 22 W/cm 2 : with 
(solid) or without (dashed) accounting for the QED effects. 
Here hu) c0 » 1.1 MeV for A = 0.8/xm. 



Radiation back-reaction. While emitting a photon, 
an electron also acquires 4-momentum from the external 
field, equal to dp F = ^^if^^ k ~ kX c (k' -pi)/{k-pi) 
(see Eg. (fTOl) ). Usually this is small compared to dp ra d- 
However, the account for the interaction with the field 
ensures that the total effect of emission on the electron 
not to break the entity (p/ ■ Pf) = 1. The choices of 
near-unity correction coefficients in dpF are somewhat 
different in the cases x < 1 and » 1. For moderate 
values of x the radiation force, (dpF — ^Prad) / 'dt, may be 
found by integrating both dpp and <ip ra d over dk': 



d{Pf ~ Pi) 
dr 



, (Pi ■ Pi 
' {k-pi] 



Iqed 
m e c z 



(15) 



where the choice of the coefficient in dpp, first, ensures 
that the radiation force maintains the abovementioncd 
entity (since (pi ■ d(p / — pi) / dr) = 0), and, second, makes 
Eq. (|15[) applicable with dimensional momenta as well. 

We already mentioned in [3] , that QED is not compat- 
ible with the traditional approach to the radiation force 



in classical electrodynamics and suggested an alternative 
equation of motion for a radiating electron: 



dp 1 
~d~r~ 



= n 



ik n 



Pk 



'QED j 
9 ¥ 



TO- 



QED olfe 



(16) 



where f2 u = eF 1 ^ /(m e c), F 1 ^ is the field tensor and 
t = 2e 2 /(3m e c 3 ). In the ID plane wave ToQ tk flkip l = 
k l (p ■ p)I c \/(m e c 2 (k ■ p)), so that the radiation force in 
Eq. lfTB"]) is the same as its QED formulation in Eq. (fl~5|) . 
This proves that the earlier derived Eq.( [T6|) has a wide 
range of applicability including an electron quasi-classical 
motion in QED strong fields. The way to solve Eq. (fj"6]) 
and integrate the emission is described in [3]. 

In Fig|4] we show the numerical result for an electron 
interacting with a laser pulse. We see that the QED ef- 
fects essentially modify the radiation spectrum even with 
laser intensities which are already achieved. 

We conclude that in a wide range of applications, 
including the case of very strong laser fields with essen- 
tial QED effects, the electron motion may be successfully 
described within the radiation force approximation. The 
necessary corrections in the radiation force and the emis- 
sion spectra to account for the QED effects are parame- 
terized by the sole parameter, I c \. 

The future application to QED Monte- 
Carlo simulations may be based on the to- 
tal probability of emission per interval of proper 
time: W = Ar(/QED/(TO e c 2 ))d; _1 , where uj^ 1 = 
J Q{ r OiX)/(x r o)drQ. The expression of the only scalar 
to gamble, Xclu'/£c, in terms of a random number, 
< R < 1, is given by an integral equation as follows: 



% c u' /(Sex) 



Q/(rox)dr = R < Xc^'/Sc > 1 



This method will be described in a forthcoming publi- 
cation in detail, including the pair production (see [l2j 
regarding the latter effect). 
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